function sse = FitExp(p, xdata,ydata )
        FittedCurve = exp_decay(p,xdata);      
        ErrorVector = FittedCurve' - ydata;
        sse = sum(ErrorVector .^ 2);
end